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Abstract 

The NASA F-106 has acquired considerable data 
on the rates of change of electromagnetic parame- 
ters on the aircraft surface during over 700 direct 
lightning strikes while penetrating thunderstorms at 
altitudes from 15 000 to 40000 ft (4570 to 12 190 m). 
These in situ measurements have provided the ba- 
sis for the first statistical quantification of the light- 
ning electromagnetic threat to aircraft appropriate 
for determining indirect lightning effects on aircraft. 
The data are presently being used in updating previ- 
ous lightning criteria and standards developed over 
the years from ground-based measurements. The 
proposed lightning standards will be the first which 
reflect actual aircraft responses measured at flight 
altitudes. Nonparametric maximum likelihood es- 
timates of the distribution of the peak electromag- 
netic rates of change for consideration in the new 
standards are obtained based on peak recorder data 
for flights which have multiple strikes. The linear 
and nonlinear modeling techniques developed provide 
means to interpret and understand the direct-strike 
electromagnetic data acquired on the F-106. The 
reasonable results obtained with the models, com- 
pared with measured responses, provide increased 
confidence that the models may be credibly applied 
to other aircraft types and used in the prediction 
of internal coupling effects in the design of lightning 
protection for new aircraft. 

Introduction 

The advent and projected application of compos- 
ite structures with flight-critical digital electronics 
compound lightning problems in advanced civil and 
military aircraft. The application of these technolo- 
gies motivates the need to quantify the lightning- 
generated environment affecting such aircraft since 
composite structures do not provide shielding equiv- 
alent to that of metal aircraft, and digital systems 
are potentially more susceptible to “upset” by elec- 
trical transients than previous analog electronic sys- 
tems. The term upset refers to the propensity for 
digital electronic systems to malfunction as a result 
of electronic transients which cause no permanent 
damage to the electronic hardware. Upsets which sig- 
nificantly affect system functional performance can- 
not be tolerated in flight-critical applications. In 
order to achieve the full performance potential of 
advanced transports or the mission potential of ad- 
vanced military aircraft employing flight-critical dig- 
ital controls and composite structures, credible de- 
sign and assessment techniques are required which 
optimally address the design of protection against 
lightning damage and upset. 


The in-flight direct-strike lightning electromag- 
netic threat criteria presently used for aircraft have 
evolved over a number of years and have generally 
been inferred from ground-based measurements of 
lightning strikes to instrumented towers. Protec- 
tion of present generation metallic transport aircraft 
against potential catastrophic effects of lightning 
through use of these criteria has been satisfactory, 
with only nuisance outages of electronic subsystems 
reported in 87 of 783 strikes (ref. 1). Application of 
composite structures, coupled with the more sensi- 
tive digital electronic subsystems, will increase the 
vulnerability of the total system function to electri- 
cal transients. The design and assessment techniques 
for future aircraft should be based on a realistic, sta- 
tistically significant definition of the lightning threat 
concerning its potential to generate interfering elec- 
trical transients. 

The NASA Langley Research Center has con- 
ducted direct-strike lightning characterization re- 
search using a specially instrumented F-106 air- 
craft since 1980. During these tests, the F-106 
has acquired direct-strike lightning data on over 700 
strikes at altitudes between 15000 and 40000 ft 
(4570 to 12 190 m). The research has emphasized the 
rates of change associated with the lightning external 
interaction with the aircraft; this emphasis was mo- 
tivated by the technical community’s interest in in- 
direct effects. The acquired electromagnetic data in- 
clude the rate of change of current to the nose boom; 
the rates of change of electric and magnetic flux den- 
sity at a number of locations on the aircraft; currents 
to the nose boom and vertical fin cap; and induced 
voltages recorded on a few internal wires. A pho- 
tograph showing the sensor locations on the F-106 
aircraft used in the Langley direct-strike lightning 
research is shown in figure 1. 

This paper includes a brief summary of the direct- 
strike data and acquisition, a description of lightning- 
aircraft interaction mathematical modeling along 
with some typical modeling results, a discussion of 
in-flight strike phenomenology, and the development 
of the statistical method used to estimate the peak 
rates of change of currents and fields. The paper 
closes with a discussion of the impact of the F-106 
data on lightning test criteria and a discussion of the 
relationship of the data to existing criteria. 

The numerical analysis presented consists of both 
linear and nonlinear finite-difference modeling. The 
analysis was performed to determine the lightning 
environment (ambient electric field, lightning cur- 
rent, channel properties, etc.) associated with each 
lightning strike to the F-106. Linear modeling makes 
several simplifying assumptions about the lightning- 
aircraft interaction to derive lightning currents. The 


nonlinear model allows the lightning-aircraft interac- 
tion to be examined in more detail, but at the ex- 
pense of far greater complexity in the model. Al- 
though the numerical models described herein deal 
with metallic aircraft only, they can be extended to 
composite structures. The effect of a composite is 
to redistribute lightning currents at later times (gen- 
erally >1 ps after the beginning of an event) and 
to allow electromagnetic energy to diffuse through 
the aircraft skin. Therefore, models of aircraft hav- 
ing composite materials must be altered to account 
for these phenomena for simulations extending be- 
yond the initial stage of a lightning event. Details of 
the models and results are presented in a subsequent 
section, 

A statistical method is proposed for estimating 
the distributions H(x) of peak rates of change of 
currents and fields on the basis of data provided by 
the F-106 peak recording instruments. Peak recorder 
data consist of pairs for i = 1,2, ...,n, where 

Wi is the largest of the individual peak readings gen- 
erated by z t strikes during the ith flight. Since the 
individually generated peak readings are not observ- 
able (because there is generally more than one strike 
on a given flight), inference concerning the form of 
H(x) must be made on the basis of a model which 
links (z t , Wj) (i — 1, 2, ..., n) to H(x). The classic way 
to deal with maxima of random variables is to con- 
sider the three possible limit laws of extreme-value 
theory (refs. 2 and 3), attempt to determine which of 
the laws is applicable to the data, and then perform 
estimation in a parametric setting. This method, 
however, requires that the number of observations 
contributing to the maxima be quite large and in- 
troduces the difficulty of discriminating between the 
three possible models. Herein a nonparametric ap- 
proach is developed in which a likelihood function 
is derived from the assumption that the individually 
generated peak rates of change are statistically in- 
dependent and have a common distribution. This is 
the first known attempt to estimate a general form 
of underlying distribution based on the observation 
of only the maxima of random variables. A compu- 
tationally feasible form of the estimating equations 
is obtained and asymptotic properties (n — > oo) of 
the maximum likelihood estimators are established. 
Details concerning the model, the method, and its 
application to the peak recorder data are presented 
in a later section. 

The goal of this research has been to study the 
lightning threat by establishing a credible lightning- 
aircraft interaction data base, obtained in situ during 
direct-strike flight tests, and by developing models 
of lightning-aircraft interaction based on the actual 
strike data obtained. This research has provided 


the first statistically significant data base on the 
peak rates of change of currents and fields on an 
aircraft when struck by lightning. The peak rates 
of change of currents on aircraft struck by light- 
ning are about two times those of previously ac- 
cepted airworthiness criteria. Since there are at 
present no criteria on the rate of change of electric 
flux density, the new data can be used as the ba- 
sis for new criteria on the electric characteristics of 
lightning-aircraft electrodynamics. The findings are 
at present being included in new criteria concerning 
protection of aircraft electrical and electronic systems 
against lightning. The new lightning-aircraft interac- 
tion modeling techniques, developed to understand 
and interpret the direct-strike measurements, can 
provide the basis for prediction of lightning-aircraft 
electrodynamics of generic aircraft types. 

Symbols 

aQ, ax, grouping limits for the observed 

largest peak rates (k -f 1 is the 
number of grouping intervals) 

matrix with elements defined by 
B rs (0) = dlps/dOr (r, 3 = 1, 2, k) 

time derivative of magnetic flux 
density 

carbon fiber composite 

diagonal matrix with the reciprocal 
of the diagonal elements being the 
variances of the limiting distribu- 
tion of n 1 / 2 ^ — t/) 0 ) 

time derivative of electric 
displacement 

thickness 

electric intensity 

mathematical expectation 

air breakdown field strength 

electric field magnitude 

electric field radial component 

relative electric field intensity 

surface electric field (tangential to 
surface) 

joint density of conditional 
on Zi = z 

avalanche rate 

magnetic field 


B 

B-dot 

CFC 

D 


D-dot 

d 

E 

E() 

E b 

E 'mag 

E r 

E rel 

E 9 

/*MM) 

G 

H 


2 


H(x) 

cumulative distribution function of 

Vi 

transmission line voltage 

individually generated peak rates of 

W: 

random variable which represents 


change of currents and fields 


the largest peak rate generated 

H(x) 

maximum likelihood estimator of 


during the zth flight 


H(x) 

Zr 

random variable which represents 

h(x) 

I 

I- dot 

h 

Is 

It 

density function of H(x ) 
electric current 

time derivative of electric current 
transmission line current 
current source 
vertical fin cap current 

Zt 
z 0 

Zi, z 2 

0 

the number of strikes to the F-106 
during the ith. flight 

transfer impedance 

characteristic impedance 

termination impedances 

electron attachment rate 

electron-ion recombination 

J 

matrix (Jrs(VO) °f second 
derivatives of the negative of the 

coefficient 

Js 

log likelihood function 
surface current density 

6 

negative-positive ion recombination 
coefficient 

k 

propagation constant 

to 

permittivity of free space 

L 

inductance 

V 

intrinsic impedance 

n 

number of flights of the F-106 

9 

parameters vector (0i, 02>-*-5 fyc+l) 


airplane 

0 

maximum likelihood estimator of 0 

n+ 

positive ion density 


— H{aj) — = 1, 2, k) 

n_ 

negative ion density 

A" 1 

covariance matrix of the limiting 

Tl e 

electron density 


distribution of n l / 2 (0 — 0 O ) 

P 

percent water vapor 

A 

electric charge per unit length 

Pi(z) 

probability that Z{ — z 

Me 

electron mobility 

Q 

ambient ionization rate 

Pi 

ion mobility 

Qm 

q 

maximum charge an aircraft may 
hold 

electric charge 

Ms 

Mr 

a 

mobility of species s 
relative air density 
electrical conductivity 

R 

resistance; response function 


n 

Rm 

measured response function 


= fi fi K'k. 0) 

i=l 

r 

radial distance coordinate 


parameters vector (^i, ^ k ) 

r trig 

radial distance at which aircraft 
lightning triggering occurs 


defined by 9\ + O 2 + ••• + 0j = 
exp ( - S ^/ ) for j — 1, 2 , k 

T 

transfer function 


V 1=3 ) 

v* 

. vector (V u , V 2i , of indi- 

$ 

maximum likelihood estimator of ^ 


cator variables defined by Vj{ — 1 
if a>j_i < W{ < dj and Vji = 0 
otherwise (j = 1, 2 , k + 1) 

0O 

true value of ip 

v s 

OJ 

Capital 

angular frequency 

letters refer to random variables and 

voltage source 

lower-case letters refer to the realizations of these 

Vs 

drift velocity of species s 

variables. 
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Direct-Strike Data Summary 

The measurements made on the F-106 are elec- 
tromagnetic parameters at the aircraft surface and 
induced voltages on a few internal wires. Several 
electromagnetic parameters have been studied: cur- 
rent and rate of change of current in the nose boom; 
current in the vertical fin cap; rate of change of elec- 
tric flux density under the forward fuselage, at the 
base of the vertical fin, and under the outboard sec- 
tion of each wing; and rate of change of magnetic 
flux density on opposing sides of the aft fuselage and 
under each wing at about midspan. The lightning 
measurements and design of the F-106 experiment 
are described in detail in references 4 and 5, and the 
instrumentation is described in reference 6. The sen- 
sors, which are generally based on designs developed 
for nuclear electromagnetic pulse measurements, are 
described in reference 7. 

The data are recorded in a shielded, self-contained 
instrumentation package which is mounted in the air- 
craft weapons bay. Power for the instrumentation is 
obtained from a motor generator set which decou- 
ples any lightning-induced transients in the aircraft 
power system to guard against spurious instrumen- 
tation system response. The electromagnetic sen- 
sors are electrically connected to the instrumentation 
package with foam-filled heliax cable. Control and 
diagnostics for the instrumentation are accomplished 
through the use of fiber-optic data links. 

The direct-strike lightning process may last about 
1 s and consists of a complex interaction of extremely 
fast electromagnetic pulses with the aircraft struc- 
ture (refs. 8 to 10). Three basic recording tech- 
niques have evolved during the research as appro- 
priate and complementary for investigating in-flight 
direct-strike lightning. The instruments used in these 
techniques are continuous analog recorders, digital 
transient recorders, and peak recorders. The con- 
tinuous analog recorders yield temporal information 
on the overall lightning process; they do not, how- 
ever, have sufficient bandwidth to record the fast 
pulses with suitable fidelity. The digital transient 
recorders have a sufficiently wide bandwidth but can 
record only during a small interval of the lightning 
event, and thus they yield information only on typical 
pulses which exceed their trigger threshold. The peak 
recorders supplement both the analog and digital 
transient recorders by obtaining information on the 
maximum parameter value attained during a strike. 

The continuous analog recorders have a nomi- 
nal bandwidth of 400 Hz to 100 kHz and, for cur- 
rent sensors with dc response, have been frequency- 
division multiplexed with subcarrier oscillators to 
obtain dc response to record the continuing current 


during strikes. The digital transient recorders, which 
have a Nyquist bandwidth of dc to 100 MHz when 
operating at the maximum selectable rate of one data 
sample every 5 ns, are central to the data acquisition 
system for obtaining time domain information appro- 
priate for development of lightning interaction mod- 
els. These recorders have 12 channels with 65 536 
8-bit words each, yielding a data window of 327 jx s 
at the maximum sample rate and have replaced the 
2-channel, 131072 6-bit word units used during the 
first several years of the research (ref. 11). The peak 
recorders sort the maximum voltage attained by the 
sensor connected to its input during a flight and 
are specified for operation over a bandwidth equiv- 
alent to half sine wave pulses with baseline widths 
between 5 ns and 10 jx s. 

Approximately 2500 individual time domain wave- 
forms have been obtained from the various sensors 
during the strikes, and 130 peak recorder readings 
have been obtained during about 400 of the strikes. 
The maximum value recorded to date on the peak 
recorder monitoring the rate of change of the nose 
boom current is 380 kA//zs; the largest rate of change 
of electric flux density under the forward fuselage 
has exceeded the full-scale range of 97 A/m 2 ; and 
the largest peak vertical fin cap current recorded 
is 54 kA. 

Figures 2 to 6 show examples of the types of time 
domain waveforms that have been recorded. These 
waveforms illustrate the overall lightning process, the 
triggering time of the transient recorder, and some 
of the wide-bandwidth transient recorder data. Fig- 
ure 2 shows the current for three different strikes 
flowing in the tip of the vertical fin cap. The strikes 
were recorded from a current sensor with dc response 
on a frequency-division multiplexed analog record- 
ing system, which also had a dc response. Note the 
unipolar dc continuing current on the order of 80 A 
lasting 0.3 s or so with a number of pulses superim- 
posed on it. As mentioned previously, the recording 
channel upper frequency response limit of 400 Hz 
does not allow sufficient fidelity to determine the 
peak values of the individual pulses, and thus the 
channel provides information useful only for study- 
ing the overall character and temporal nature of the 
lightning currents. 

Figure 3 shows the temporal character of a typical 
lightning strike obtained from a number of different 
sensors recorded with a continuous analog recorder 
with a bandwidth of 400 Hz to 100 kHz. (The 
fin cap current was also recorded with a dc-to-400- 
Hz bandwidth channel.) This figure shows pulses 
from a light sensor which are concurrent with pulses 
from other sensors measuring the lightning external 
interaction with the aircraft. (The pulsing nature of 
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the light from lightning attached to the aircraft has 
also been noted in data from cockpit-mounted video 
cameras.) 

Figure 4 shows details of the fin cap current of 
figure 3 recorded at a 40- ns sample interval with a 
digital transient recorder at the time of the transient 
recorder trigger. The trigger time is shown as the 
discrete event at the bottom trace of figure 3. The 
current in figure 4 has a high repetition rate, with 
about 50 //s between the two 12-kA pulses. 

Figures 5 and 6 are typical simultaneous digital 
transient recorder waveforms (recorded at a 10-ns 
sample interval) of rates of change of electric and 
magnetic flux density from sensors located under the 
forward fuselage and on the side of the aft fuselage. 

A multitude of time domain waveforms have been 
published which include all the transients recorded 
through the 1984 tests for 627 strikes (refs. 12 to 16). 
Photographs of lightning attachments taken from 
onboard motion picture cameras during 156 strikes 
in 1982 are in reference 17. Strike conditions, at- 
tachment photographs, and lightning patterns on 
the aircraft structure are discussed in reference 18. 
An overview of the acquired data including selected 
waveforms is given in reference 19. 

Modeling and Phenomenology of Lightning 
Strikes to Aircraft 

Lightning-Aircraft Interaction Modeling 

In order to understand the lightning data col- 
lected from the F-106 thunderstorm penetrations and 
to extend that understanding to the lightning en- 
vironment, it is necessary to model the lightning- 
aircraft interaction. Several distinct models have 
been developed in the course of this investigation. 
These include both nonlinear models for analyzing 
the physics of a lightning event and linear models 
(ref. 19) which can analyze the interaction in a simpli- 
fied, or engineering, sense. All these models are based 
on the finite-difference methodology first developed 
by Yee (ref. 20). 

In this section, the methodology of the finite- 
difference modeling is discussed briefly, followed by 
specific descriptions of the linear and nonlinear mod- 
eling techniques and selected results from these mod- 
els. Comparisons with measured flight data are 
included. 

Three-Dimensional Finite-Dif ference 

Modeling 

The main analytical tool used in the analy- 
sis of the lightning-aircraft interaction is the computer 


code T3DFD, which stands for time domain 
3-dimensional finite-difference code (ref. 21). The 
code solves Maxwell’s equations in three dimen- 
sions and is capable of modeling complex geometries, 
space- and time-varying permittivity, permeability, 
and conductivity. It has been enhanced with an air 
chemistry model to allow the calculation of air con- 
ductivities which occur in the vicinity of large electric 
fields. 

The particular problem space used in the F-106 
program is a Cartesian mesh enclosing a space ap- 
proximately twice the size of the aircraft itself. The 
spatial resolution is 1 m in the direction along the 
fuselage and 0.5 m in the wing-wing and vertical 
directions. The temporal resolution of the model 
is 1 ns. The nominal frequency resolution of the 
mesh, if we assume a minimum of five cells per wave- 
length, is therefore 60 MHz. The F-106 is placed 
within the problem space with all electric fields tan- 
gential to the surface set at zero at all times. Hence, 
the aircraft is assumed to be perfectly conductive, 
with no significant apertures which would alter the 
external response. This is a good approximation 
except for the immediate vicinity of the cockpit. 

The F-106 as it appears in the finite-difference 
code is shown in figure 7. The coordinate system is 
also indicated, as are the locations of the external 
sensors which are modeled. One thing to notice 
about figure 7 is the blocky nature of the F-106 
model. The large-scale structure of the aircraft is 
well resolved, but it is clear that details, such as 
the nose boom, are not resolved. This can cause 
difficulty when field distributions around such points 
are desired. In most cases, however, the model as 
shown is adequate to predict the response of the 
sensor system to a given lightning event. This is 
because the sensors have been placed away from 
those portions of the plane that are not well resolved 
by the finite-difference code. 

Nonlinear Modeling 

The development of a lightning strike to an air- 
craft certainly involves nonlinear processes. In par- 
ticular, the initiation of a triggered lightning event 
by an aircraft immersed in a high electric field must 
include the development of a corona and associated 
conductivity around the aircraft, which is a func- 
tion of the surrounding electric field. In order to 
correctly predict nonlinear aircraft responses to trig- 
gered lightning, any nonlinear model must be able to 
adequately model the corona growth. 

The model described herein for an electrical 
corona solves for the air conductivity. It does this 
by calculating the densities of positive ions, negative 
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ions, and electrons as a function of space and time 
through the use of detailed balancing. Physical pro- 
cesses included are electron avalanching, electron at- 
tachment to neutral molecules to form negative ions, 
electron-positive ion recombination, and negative- 
positive ion recombination. 

The air conductivity for the model is a nonlinear 
function of the electric field magnitude and can be 
written as 


o = q{n e He + {n++n-)ni) ( 1 ) 


Here, q is the electric charge, n e is the free electron 
density, n_ is the negative ion density, n+ is the 
positive ion density, fi e is electron mobility, and fii is 
the ion mobility. 

The electron and ion densities are calculated from 
a rate equation for each species: 


+ V * (™eV e ) + + a e — G] n e — ^ (2) 


dt 


+ V • (n_v_) + 6n+n- — a e n e 


( 3 ) 


0 ^ ~t~ V • ( n + v +) + (3n e n+ + 6n-n-\. — - T Gn e (4) 


where is the ambient ionization rate, a e is the 

electron attachment rate, / 3 is the electron-ion recom- 
bination coefficient, and 8 is the negative-positive ion 
recombination coefficient. The velocity v«* of each 
species s is determined from the mobility and the 
electric intensity E: v s = /i s E. The actual rates used 
in the model are analytical fits to experimental data 
and are given in appendix A. 

The closure of the rate equations for each species 
with the mobility equation for the velocities is a 
step which ultimately limits the model to low- 
temperature phenomena. Ideally, the set of equations 
should be supplemented by momentum and energy 
conservation equations and closed by an appropriate 
equation of state. This would allow the model to fol- 
low the development of a lightning event for a much 
longer time. 

The integration of the nonlinear corona model 
with T3DFD is accomplished through the use of the 
crE term in Maxwell’s equations (refs. 22 to 25). 
The T3DFD code provides an electric field magni- 
tude at each point in space; from the magnitude 
the air chemistry model calculates a conductivity 
which is then fed back into Maxwell’s equations 
through the crE term. The air chemistry equations 
are also solved with a finite-difference method, al- 
though somewhat differently than Maxwell’s equa- 
tions are solved. Maxwell’s equations are central 


differenced in time and space for maximum stabil- 
ity. The air chemistry equations, because of their 
nonlinear character, cannot be efficiently central dif- 
ferenced and therefore one-sided differencing is used. 
Numerical tests with varying space and time steps 
have shown that for these equations, one-sided dif- 
ferencing is sufficiently accurate and does not affect 
overall stability of the model. 

One of the uses of the calculated responses from 
the parameter study is for comparison with measured 
responses. If one can find an element of the parame- 
ter study that matches the measured data for a par- 
ticular lightning event, then there is reasonable con- 
fidence that the ambient environment responsible for 
the event was similar to that of the parameter study 
element. Because of the nonlinear character of the 
parameter study, no scaling of calculated response 
amplitudes can be done, which implies that compar- 
isons can be done in a general sense only, without 
detailed matching of amplitudes. The comparisons 
presented in this section were chosen mostly on the 
basis of waveform comparison without regard to am- 
plitude, although in some cases the amplitudes also 
match quite closely. Examples of comparisons are 
shown in figures 8 to 11. 

Of particular interest in the comparisons is the 
double pulse character shown in figures 8, 9, and 11. 
This double pulse corresponds to an air breakdown at 
two different places on the F-106 at separate times. 
The first breakdown raised the electric field at the 
second location enough to cause a breakdown there 
also. This type of comparison is strong evidence that 
the triggered lightning events on the F-106 are well 
modeled by the nonlinear air breakdown model. 

Linear Modeling 

Nonlinear lightning modeling, although more 
physically complete than linear methods, is costly 
and time-consuming. It is therefore quite useful to 
develop linear modeling techniques which allow one 
to investigate certain aspects of the lightning event 
without resorting to the complexity of a nonlinear 
model. One such technique is the use of transfer func- 
tions to determine the lightning current which caused 
a given set of measured electromagnetic responses on 
the F-106 aircraft. The transfer function is a func- 
tional relationship in the frequency domain between 
a source function and a response function. The use 
of a transfer function requires that the system under 
consideration be linear. This requirement is satisfied 
by a linear finite-difference code, but of course it is 
not satisfied by a real lightning-aircraft system. The 
justification for using the transfer function technique 
is that the nonlinearity in the real system is confined 
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mostly to the lightning channel itself. The electro- 
magnetic responses on the aircraft are often approxi- 
mately linear functions of the lightning current which 
flows onto the aircraft at the lightning attachment 
point. That is, although the formation of the light- 
ning channel, its evolution, and the lightning cur- 
rent are complicated nonlinear functions of geometry 
and initial conditions, the aircraft responses usually 
depend, in a linear fashion, only on the current at 
the attachment point. 

There are a number of assumptions which must be 
made when the transfer function technique is used. 
These assumptions are discussed individually below. 

1. Lightning attachment points must be known. 
These attachment points constitute part of the initial 
conditions and as such are necessary to define the 
problem. In addition, the attachment points cannot 
change with time because time-dependent geometry 
is not presently allowed in the code. Also, if there are 
multiple channel attachments, only one of these can 
act as a source while the others just drain charge from 
the aircraft. Although having more than one source 
does not violate linear constraints, the problem no 
longer has a unique solution if more than one source 
is involved. 

2. Relative formation times of multiple channels 
must be known. Because it is likely that exit channels 
for lightning currents form later than entry channels, 
it is necessary to know their formation times. In 
a sense, these channels which appear during the 
course of a problem constitute boundary conditions 
which change with time. This does not violate the 
linearity requirement as long as the formation times 
are fixed and do not vary with the time evolution of 
the problem. 

3. Lightning channel geometry should be known. 
Although less important than the first two require- 
ments, it is desirable to know the orientation of the 
lightning channel with respect to the aircraft. This 
geometry also constitutes boundary conditions which 
can affect aircraft responses to some extent. The re- 
sponses are affected because electromagnetic radia- 
tion from current in the channel also produces some 
response on the aircraft, in addition to the current 
which flows onto the aircraft. The contribution of 
the current is much larger, however, so the channel 
radiation can be considered as a perturbation. 

4. Lightning channel impedance should be known. 
Again, this is less important than the first two re- 
quirements, but it can affect the linear relationship 
between source and response somewhat. The reason 
for this is that the back effect of the aircraft on the 
lightning current is different depending on the chan- 
nel impedance. For example, a channel with infinite 
impedance is completely unaffected by the presence 


of the aircraft. In this case, whatever current is flow- 
ing in the channel is the current which is transferred 
onto the aircraft, and reflections and responses on 
the aircraft behave as if the channel were nonexis- 
tent. For any other impedance, the response of the 
aircraft can affect the current in the channel, espe- 
cially for lower impedances. The lower impedance 
slightly alters the electromagnetic responses of the 
aircraft. Hence, it is desirable to know the channel 
impedance as a boundary condition of the problem. 
It is possible to handle a time-varying impedance, as 
long as the variation is specified in advance and does 
not depend on the evolution of the problem. 

The transfer function technique can be used for 
either triggered or natural lightning. The distinction 
between the two types of lightning is that triggered 
lightning begins at the aircraft and moves away, 
whereas natural lightning begins away from the air- 
craft and moves toward it. Typical geometries for 
both of these cases are shown in figure 12. The 
case for natural lightning is shown in figure 12(a), 
A current channel having a specified impedance and 
velocity of propagation is attached to the nose of the 
F-106. The impedance and velocity of propagation 
are determined from the inductance per unit length 
and the capacitance per unit length of the channel, 
which are determined from the physical diameter of 
the channel. Therefore, in the linear model, all the 
electrical properties of the channel are determined 
from the size of the channel. The attachment to the 
nose depicted in figure 12 is simply an example of the 
point at which many lightning strikes attach to the 
F-106. The model allows attachment at any point on 
the aircraft and allows multiple channels as long as 
there is only one current source in the problem. The 
difference between figures 12(a) and 12(b) is in the 
location of the current source. For natural lightning, 
the current source propagates toward the aircraft and 
is located at the edge of the problem space, in order 
to model initiation and driving forces for a natural 
lightning strike removed from the aircraft. The cur- 
rent source in figure 12(b) is located near the point 
where the channel attaches to the F-106 in order to 
model the phenomenon of initiation at the surface of 
the aircraft with propagation outward. 

The mathematics of the model are handled iden- 
tically in the two cases. A transfer function is de- 
termined from the Fourier transforms of the current 
source and the response waveforms: 

TV) = *V)//M (5) 

where T(uj) is the transfer function in the frequency 
domain, R( u) is the Fourier transform of the cal- 
culated response waveform, and I(oj) is the Fourier 
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transform of the current source waveform. Since the 
model is linear, T(u/) is source independent; that is, 
changing I changes R in such a way as to keep T the 
same. 

To determine the current source necessary to pro- 
duce a given measured response, the measured re- 
sponse waveform must be Fourier transformed. Then 
the transfer function for the particular sensor is used 
to calculate the Fourier transform of the current 
source necessary to produce the measured response 
with equation (6). 

/(«) = iZmM/TV) (6) 

The current is then transformed into the time domain 
where it can be used in the linear model as a check 
to ensure the measured response is reproduced. 

The transfer function technique is most useful in 
the events for which multiple simultaneous sensor re- 
sponses have been recorded. Several current sources 
can then be derived, one for each of the multiple re- 
sponses. In principle, if the model geometry chosen 
is correct, all the current sources should be the same. 

An example of results from the linear model is 
presented in figure 13. The model used is of a 
triggered lightning strike to the nose of the F-106 
with an exit channel also forming at the tip of 
the vertical stabilizer. The transfer function used 
to derive the source current is that appropriate to 
the B - dot longitudinal sensor. Agreement between 
predicted and measured results for the other sensors 
is quite good, particularly for the left wing sensors. 

Phenomenology of Lightning Strikes to Aircraft 

The phenomenology of an aircraft encountering 
a lightning event is complex and involves several 
parameters. These parameters include aircraft size 
and shape, altitude, static charge, ambient electric 
fields in which the aircraft is immersed, location of 
charge centers, and particle environment. 

In order to discuss this phenomenology, it is con- 
venient to categorize aircraft-lightning interaction 
events into three types: natural lightning, triggered 
lightning, and hybrid lightning. It is our intent in this 
section to define these concepts and to describe the 
basic phenomenologies associated with each category. 

Natural Lightning 

Natural lightning occurs without the presence of 
the aircraft. The presence of the aircraft is only 
incidental to the occurrence of natural lightning and 
is not contributory. An aircraft is subjected to this 
environment when, by chance, it is in the path of a 
naturally occurring lightning strike. 


This phenomenon is further described with the 
help of figure 14. Figure 14(a) shows a stepped leader 
propagating from a cloud toward the Earth, and the 
aircraft happens to be in its path. Because the ap- 
proaching leader induces large electric fields on the 
aircraft, arcs develop on the aircraft extremities, and 
these arcs propagate toward the leader and make 
contact, as shown in figure 14(b). The leader may 
then continue from the aircraft toward the ground, 
as in figure 14(c), or another branch of the channel 
may propagate toward the ground, as in figure 14(d). 
In either case, the important observation is that the 
lightning channel would have existed without the air- 
craft, and the aircraft’s presence did not significantly 
influence the path of the channel. 

The same scenario would apply to a natural in- 
tracloud event. It is believed that the F-106 has 
encountered very few natural lightning events. 

Aircraft-Triggered Lightning 

The focus of this paper is on aircraft-triggered 
lightning because most of the strikes to the F-106 are 
triggered events. This assessment is based largely on 
the observations in reference 26 describing the anal- 
ysis of 49 radar returns of the F-106 being struck by 
lightning. In every case, the data in reference 26 in- 
dicate the channel begins at the aircraft and travels 
away from it. (It was later found that most, but not 
all, strikes were triggered (ref. 27).) This evidence 
prompted the development and application of non- 
linear plasma physics air breakdown modeling of the 
aircraft triggering previously described which yields 
results in good agreement with measured data. 

Aircraft-triggered lightning occurs because of the 
presence of an aircraft and would not otherwise ex- 
ist. It occurs because the aircraft is able to suffi- 
ciently enhance the local electric field to cause air 
breakdown. The largest enhancements occur at the 
sharp points or edges, particularly if those points and 
edges are oriented in the direction of the ambient 
field. An aircraft in flight, such as the F-106, has 
many sharp points (e.g., the nose) and sharp edges 
(e.g., the wings and tail) around which the fields will 
be enhanced if the aircraft is immersed in a static 
electric field. These locally enhanced fields (within 
0.5 m or so from the aircraft) are likely to be at 
least 10 times larger than the ambient field. Hence, 
it is considerably more likely that initial air break- 
down and formation of a lightning channel will oc- 
cur in the presence of the F-106 than in its absence. 
The requirements for triggered lightning are simply 
stated as follows. The static electric field which the 
aircraft experiences must be large enough and ori- 
ented properly so that the locally enhanced fields 
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at the extremes of the aircraft exceed the local air 
breakdown value. The ambient field must also be 
large enough to form and propagate a channel once a 
breakdown region has been established. A net charge 
on the aircraft also contributes to the local fields and 
can either suppress or enhance the breakdown, which 
may significantly change the requirements for the oc- 
currence of triggered lightning. Also, the net charge 
can change the location of initial air breakdown by 
altering the local field distribution around the air- 
craft. Finally, one may conclude that triggered light- 
ning will occur, if the proper conditions are satis- 
fied. Essentially, the local field somewhere on the 
aircraft must exceed the air breakdown values of ap- 
proximately 1.5 MV/m at 20000 ft and 3.0 MV/m 
at sea level. In addition, triggered lightning will oc- 
cur at the lowest possible field level. Hence, all the 
initial electromagnetic transients should be similar in 
amplitude depending on the location on the aircraft 
at which the initial breakdown occurs. A more de- 
tailed look at the environment required for triggered 
lightning is presented in the next section. 

Triggered lightning environment . In order to 
discuss the* triggered lightning environment, a typical 
thunderstorm is used as an example (ref. 28). The 
static electrical characteristics of this thunderstorm 
are a 40-C charge centered at 10 km above ground, 
a —40-C charge centered at 5 km, and a 10-C charge 
centered at 2 km, as shown in figure 15. The model 
of figure 15 is simplified slightly if the charges are 
assumed to be point charges located in a vertical line. 
For these simplifications, the static electric field can 
be analytically calculated as a function of space. 

The calculation assumes that the Earth is a per- 
fect ground plane, so it can be replaced by image 
charges. The results of the calculations are shown in 
figures 16(a) to 16(c). These figures show contours of 
a constant electric field as a function of space. The 
vertical scale is altitude and the horizontal scale is 
radial distance from the (assumed) vertical line of 
charges. Figure 16(a) shows contours of the radial 
component of the field, figure 16(b) shows the ver- 
tical component of the field, and figure 16(c) shows 
the total field magnitude. It should be remembered 
that the contours shown on the figures were calcu- 
lated with the assumption of point charges, so the 
fields shown are actually upper limits on the actual 
field. If one of the contours is inside an extended 
charge cloud, the actual field is likely to be somewhat 
smaller at that point than the contour indicates. 

Figure 16(a) is the most meaningful for the pur- 
poses of investigating triggered lightning on the 
F-106. If the aircraft is always in essentially level 


flight, the horizontal component of the field will be 
enhanced the most. 

The triggered lightning response of the F-106 has 
been successfully modeled by the nonlinear approach 
previously described. Appendix A gives the air 
chemistry coefficient formulas used in the nonlinear 
model. Inherent to this approach is the calculation 
of the static field enhancements about the aircraft as 
shown in appendix B. Electrostatic field calculations 
done on the F-106 show that fields along the axis of 
the aircraft fuselage are enhanced by about a factor 
of 10 at the nose. Wing-wing ambient fields are 
enhanced by about a factor of seven at each wing 
tip. Vertical ambient fields are enhanced by about a 
factor of three at the tip of the vertical stabilizer. 
Hence, for level flight, vertical ambient fields are 
much less important than radial fields, unless the 
vertical ambient field is several times larger than the 
radial field. 

Another factor that must be considered is that 
the minimum air breakdown field is larger at lower 
altitudes because of the increased air density there. 
A field which causes triggered lightning at a high 
altitude (^10 km) may not do so at lower altitudes 
because of the greater air density. The dashed lines in 
figures 16(a) and 16(b) illustrate this effect. They are 
the locus of points in space at which an electric field 
of breakdown strength can be reached on the aircraft 
with proper orientation. We assume, of course, that 
the net aircraft charge is zero. For a nonzero charge, 
the dashed lines could extend farther away from the 
charge centers. They also indicate the regions inside 
which triggered lightning can occur for the F-106 in 
level flight. Note that the volumes are much larger 
at high altitude, extending to almost 2 km from the 
charge center for the radial field case. At the very 
low altitudes (^2 km) the trigger region is only a 
few hundred meters across. In fact, because the 
charges are really extended and not point sources, 
it is possible that at these low altitudes no triggering 
region exists at all. 

The response of aircraft to triggered lightning . In 

order to interpret direct-strike data acquired with the 
F-106, it is of some interest to investigate in a simple 
way the expected response of the aircraft. That is, 
given that triggered lightning occurs, what type of 
measured surface 5- dot and 5-dot records can be 
expected? 

To begin, consider the situation shown in fig- 
ure 17. This situation corresponds to the F-106 flying 
directly toward a positive charge center. The air- 
craft is polarized oppositely from the ambient field, 
resulting in a field at the 5-dot sensor point (under- 
neath the nose), which points toward the fuselage. 
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Presumably this field grows slowly over several sec- 
onds as the aircraft flies into a slowly increasing am- 
bient field. This slow growth is far below the trigger 
level for the derivative sensors on the aircraft. For 
this geometry, the largest enhanced fields are at the 
nose of the aircraft, and eventually the enhanced field 
there becomes large enough to cause a corona. At 
that time electrons flow off the nose causing the field 
at the B- dot sensor location to become less negative. 
This happens very quickly and probably triggers the 
onboard instrumentation. As the air breaks down 
at the front of the F-106, the electric field increases 
at the aft end, possibly causing another breakdown 
there. That event is dependent on the characteris- 
tics of the breakdown at the nose. The corona at 
the nose effectively increases the length of the F-106 
and therefore enhances the local field at the rear of 
the aircraft. Therefore, a large corona at the nose is 
more likely to cause a breakdown behind the aircraft. 

Figure 17 also shows the expected behavior of the 
fields and their time derivatives at the B-dot forward 
and B-dot longitudinal sensor positions. The plots 
are simplified and ignore any behavior caused by the 
resonances of the aircraft. Actual records could be 
expected to look like those in figure 17 only with 
respect to general features. The pulses seen there 
could be separated by several hundred nanoseconds 
or longer, depending on the growth of the corona at 
the nose. It is significant that both single pulses and 
double pulses, as shown in figure 17, are seen in the 
measured data. 

An example from the measured data that illus- 
trates this behavior is shown in figure 18. Presented 
therein are two B-dot and B-dot events separated by 
about 170 fis. The waveforms behave similarly to the 
simplified waveforms shown in figure 17. 

The case in which the ambient field is reversed 
in direction is similar. This reversal would occur if 
the aircraft were flying away from a positive charge 
center or toward a negative charge center. If one 
assumes once again that the initial breakdown is at 
the nose and a possible second breakdown occurs at 
the aft end, the expected field behavior is inverted 
from the previous case. The actual situation is some- 
what different for the two cases, however. In the case 
shown in figure 17, a negative corona forms around 
the nose because the front of the aircraft is nega- 
tively charged. In the other case, a positive corona 
forms there. It is well-known that a negative corona 
requires a higher electric field to form than a positive 
corona, so the details of the two cases can be expected 
to differ somewhat (ref. 29). Intuitively, it seems rea- 
sonable that the smaller field required to initiate a 
positive corona would result in a slower and smaller 
measured B-dot response as the corona expands 


outward from the aircraft. The negative corona re- 
sults in a faster breakdown when the required higher 
field is finally reached. This appears to be charac- 
teristic of the measured data, in that initially posi- 
tive B-dot records (near the F-106 nose) tend to be 
larger and to have more high-frequency content than 
initially negative records. 

This very simple model illustrates the environ- 
ment for triggered lightning in the vicinity of a typ- 
ical thunderstorm. In the discussion so far, three 
complicating factors have been omitted. First, in a 
real thunderstorm, the charge centers are not points 
but exist in extended regions. They also may be 
clumped into several smaller charge centers, around 
which locally high electric fields may exist. Second, 
in most thunderstorms there is a shear in the verti- 
cal direction, so that the charge centers are not in a 
vertical line. This destroys the cylindrical symmetry. 
Third, other thunderclouds in the vicinity alter the 
field distribution around any given cloud. 

Fields required for triggering the F-106 at high 
(8230 m) and low (3050 m) altitudes have been calcu- 
lated in a parameter study for various aircraft orien- 
tations and values of net charge on the aircraft. Light- 
ning triggering conditions are determined by first 
computing the static fields on the aircraft caused by 
the net charge. The fields on the aircraft caused by 
the static ambient thunderstorm fields can be added 
through linear superposition to the fields caused by 
the net charge. The trigger level is then the ambient 
thunderstorm field magnitude, which causes the to- 
tal field on the aircraft to exceed the corona level at 
at least one location. The procedure is to place the 
F-106, with a net charge on its surface, in a uniform 
electric field through the use of linearly derived ini- 
tial conditions. The magnitudes of the field and the 
charge are chosen so that an air breakdown of the 
same order of magnitude as that seen in the measured 
data occurs (i.e., calculated and measured B-dot and 
B-dot responses on the aircraft are of similar mag- 
nitude). The location of this breakdown varies de- 
pending on the orientation of the field and the sign 
of the charge, but it generally occurs in one of three 
places: the nose, one of the wing tips, or the tip 
of the vertical stabilizer. Occasionally two locations 
are involved, as a breakdown at one location raises 
the local field at another location enough to cause a 
breakdown at the second location also. 

The results of the triggering field calculations are 
summarized in table 1. Along the upper row of 
the headings in table 1 is shown the orientation of 
the ambient electric field for each of the parameter 
study elements in that column. The orientations are 
along the principal coordinate axes of the aircraft — 
longitudinal, horizontal, and vertical. Both positive 
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and negative fields are considered, because a negative 
corona air breakdown is different than a positive 
corona air breakdown. The aircraft charge Qm is 
computed to be 1.79 mC and is the value for which a 
corona would form on the aircraft by virtue of only 
its net charge. The Q m is thus the maximum charge 
the aircraft can hold. In the leftmost column is the 
net charge which was placed on the aircraft for each 
orientation. Five charge values were used: no charge 
(0), ±0.5Q m , and ±Q m . 

The field magnitudes for a given charge and orien- 
tation are given in the columns of values. Note that, 
as expected, the field magnitudes are much less for 
the cases having a charge, because the charge pro- 
duces virtually all the local fields necessary to pro- 
duce air breakdown. A word of explanation is neces- 
sary as to why these cases require any ambient field 
at all. Because Q m was chosen to be large enough 
to cause breakdown by itself, it would seem that any 
field is superfluous. But this is not always the case; 
for example, a positive Q m charge on the aircraft in 
the absence of a field induces air breakdown at the 
nose. However, if even a small ambient field oriented 
nose to tail is present, the breakdown at the nose is 
suppressed. For that field orientation and that net 
charge, the field magnitude must be increased until 
breakdown occurs at some point other than the nose. 
An example of the calculated sensor outputs for one 
element of the parameter study is shown in figure 19. 
This includes the outputs A to D. A typical nonlinear 
current (output E) is shown in figure 20. 

There are six outputs from each element of the 
parameter study. The first is the electric field mag- 
nitude necessary to cause air breakdown for a given 
charge on the aircraft and a given orientation of the 
field (table 1). The other five outputs are responses 
of the aircraft to the initial conditions. These are 
listed below: 


Designation 

Output 

A 

D-dot forward 

B 

D-dot tail 

C 

D-dot left wing 

D 

B-dot longitudinal 

E 

Current at suspected air breakdown 
point 


The locations of outputs A to D are shown on the 
finite-difference block model of the F-106 in figure 7. 


The location of output E was chosen by estimating 
where breakdown would occur. The estimate was 
made by locating the place on the model at which 
the largest electric field existed. This was the point 
of initial corona formation, but in some cases not the 
point where the largest current flowed at later times. 
Therefore, output E was not always useful. 

In addition to the principal angles of electric 
field orientation with respect to the F-106 described 
above, the parameter study has been extended to 
intermediate angles between those shown in table 1. 
This was done for the higher flight altitude only. The 
results are similar to those already presented and are 
therefore omitted. 

The effects of aircraft size and shape on trigger 
conditions . The overall objective of the F-106 direct- 
strike research program is to determine lightning 
environments for any size or type of aircraft. By 
necessity, nonlinear analysis in this program has con- 
centrated on the F-106 itself. The thunderstorm en- 
vironment necessary for a direct strike may be dif- 
ferent, however, for different types of aircraft. That 
is, for an identical thunderstorm environment, one 
aircraft may encounter a direct strike, while a sec- 
ond differently shaped or sized aircraft may not. For 
this reason, it is desirable to look at the environment 
more generally. In this section, five different aircraft 
are placed in a thunderstorm environment (i.e., am- 
bient electric field, charge on the aircraft, particu- 
lar altitude, etc.) and the conditions necessary to 
produce a triggered lightning strike are determined. 

The five aircraft investigated are listed below: 

1. Half-size F-106 

2. Double-size F-106 

3. Normal-size F-106 with delta wings 
replaced with more conventional 
straight wings 

4. Normal-size C-130 

5. Normal-size F-106 

The half- and double-size F-106 aircraft models 
afford a size range from approximately 9 to 36 m. Di- 
mensional drawings of the scaled aircraft are shown 
in figures 21(a) to 21(e). Although the field enhance- 
ments around these two aircraft are, of course, the 
same, the effect of a given net charge on the aircraft 
is much different. Also, the responses of electromag- 
netic sensors on these aircraft are different because of 
the significantly different resonance characteristics. 

The F-106 with the straight wings is used to eval- 
uate the effect of the delta wing on field enhance- 
ments, resonant response, and charge storage prop- 
erties. These are all expected to be changed. In 
particular, charge storage could be significantly af- 
fected, because the presence of the large, flat delta 
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wings allows the F-106 to hold more charge before 
breakdown field strengths are achieved. 

The C-130 is included as a nonfighter type of 
aircraft. It has more rounded contours, particularly 
at the nose, and is representative of the shapes found 
in large passenger aircraft. It should give an idea 
of the environment necessary to generate triggered 
strikes on that type of aircraft. The normal-size 
F-106 is included for comparison purposes and for 
determination of relationships with actual flight data. 

For each aircraft, two different environments were 
used. In both environments the net charge on the 
aircraft was fixed at — 0.5Q m , where Q rn is the max- 
imum charge which can be on the airplane without 
causing immediate corona formation. It should be 
kept in mind that Q m is different for each aircraft. 
Two separate ambient field orientations were used, 
one with the field oriented from nose to tail and the 
other oriented from top to bottom. For consistency, 
all computations were done at a simulated altitude of 
8230 m (relative air density of 0.5) and with a water 
vapor percentage of zero. 

The results are summarized in table 2. The field 
values E mag represent the minimum field necessary 
to produce a triggered strike on the given aircraft 
model for the given orientation and Q m . The Q m is 
listed in millicoulombs and E m ag is in kilovolts per 
meter. 

Several noteworthy items are apparent from ta- 
ble 2. First, aircraft size is the dominant determin- 
ing factor for Q m . For aircraft with similar shapes, 
one would expect Q m to scale as the square of a typ- 
ical dimension, along with the surface area. This is 
shown to be true, as Q m varies by a factor of ap- 
proximately 16 for the factor of 4 difference between 
the half-size and the double-size F-106. Second, the 
effect of the delta wing is shown in the table. The 
magnitude of Q m for the normal-size F-106 with a 
delta wing is 1.79 mC versus a value of 2.00 mC for 
the F-106 with a straight wing. This is a somewhat 
surprising result, as one would expect that the large 
surface area of the delta wing would allow it to hold a 
large amount of charge. However, the presence of the 
delta wing changes the static electrical characteris- 
tics of the aircraft considerably. Because of the delta 
wing at the rear of the aircraft, field enhancements 
are significantly larger around the nose. In fact, a 
nose-to-tail oriented field has an enhancement factor 
of 6.3 for the aircraft with a straight wing and 7.4 for 
the aircraft with the delta wing. The same type of 
phenomenon occurs for a given net charge. That is, 
the charge arranges itself on the delta wing aircraft 
in such a way that it takes less charge to bring about 
corona formation. 


Of importance also in table 2 is the magnitude of 
the field required for a lightning strike on the C-130. 
For the nose-tail field orientation the field magnitude 
is appreciably larger than that required for the F-106 
models, which can be attributed to the bluntness of 
the nose of the C-130. In actual practice, the trigger- 
ing field may differ by even more. The reason for this 
is that the C-130 is quite accurately represented near 
the nose by the finite-difference model, so its maxi- 
mum enhancement should be correct. However, the 
sharp nose boom on the F-106 is less accurately rep- 
resented, so actual enhancements on the F-106 are 
likely to be larger than the model predicts. 

Hybrid Lightning 

In addition to natural and aircraft-triggered light- 
ning, a third type of strike to an aircraft is possible 
which will be called a “hybrid” strike. In this case, a 
natural leader channel is positioned somewhere near 
the aircraft, but not near enough for the aircraft field 
distortion to alter the channel course. The lightning 
leader leaves behind it in the channel a certain charge 
per unit length which produces, to the lowest order, 
a static field around the channel. This static field, 
at the location of the aircraft, may be large enough 
to cause a triggered streamer from the aircraft. The 
streamer would then certainly propagate to the natu- 
ral channel, becoming a branch of the natural strike. 
It is unclear whether a strike of this kind would be 
more characteristic of a natural strike or a triggered 
strike. 

A simple calculation can be done to get some idea 
of how near the aircraft a natural strike must be to 
cause a hybrid strike. For simplicity, one may assume 
that the natural channel is vertical and has a uniform 
charge density per unit length. Because only the 
orders of magnitude are desired, it is also assumed 
that the situation is quasi-steady state. That is, 
changes in the field distribution around the channel 
occur slowly and the static solution is approximately 
correct. The electric field E r around a channel charge 
per unit length A, which is effectively infinitely long, 
is radially directed and is given by 


Choosing an enhancement factor for the aircraft 
of about 10 , which is typical around the nose of the 
aircraft, one finds that triggering occurs at a distance 

= {8) 

where E 5 represents the air breakdown field strength. 
A typical value for A, calculated for a channel 
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having a charge of 5 C over a 3-km distance, is 
1.667 C/km. The value E ^ varies with altitude; at 
sea level it is about 3.0 MV/m and at 8230 m it is 
about 1.5 MV/m. Substituting these numbers into 
the above equations gives triggering distances of ap- 
proximately 100 m at sea level and 200 m at an 
8230-m altitude, resulting in a significantly greater 
distance than the distance at which the presence of 
the aircraft can directly alter the course of a natural 
lightning channel. Hence this hybrid strike is much 
more probable than a natural strike. It may even 
be that the natural strike as defined here does not 
exist, in that all strikes involving a natural lightning 
channel occur according to the hybrid scenario. 

It is possible that the hybrid process could ac- 
count for some of the variability in the measured data 
on the F-106. The reasoning behind this is the fol- 
lowing. The normal triggered lightning model, which 
has the aircraft flying into a slowly increasing ambi- 
ent electric field, requires that the triggering occur at 
the lowest possible threshold field level and results in 
aircraft responses which are very similar in amplitude 
and structure, because the only really significant ad- 
justable parameters are the orientation of the field 
with respect to the aircraft and the net charge on 
the aircraft. This smallest field level is the electric 
field strength at which air breakdown occurs for any 
fixed altitude. Above this level air breakdown will 
always occur. Because an aircraft’s motion through 
cloud fields is slow compared with electrodynamic re- 
laxation time scales and air breakdown time scales, 
air breakdown must always occur at this threshold 
field level. However, in the hybrid model the natural 
channel can appear at any distance from the aircraft 
in a relatively short time, so triggering may occur 
at field levels which are significantly above the mini- 
mum level, resulting in aircraft responses which have 
much wider variety in both amplitude and structure. 

Given the conditions under which natural and 
triggered lightning strikes occur, it appears far 
more likely that most strikes to the F-106 are trig- 
gered. Triggered lightning is certain to occur un- 
der the proper conditions, while natural lightning is 
probabilistic even under ideal conditions. 

Nonparametric Maximum Likelihood 
Estimation of the Distributions of Peak 
Rates 

The magnitudes of the voltages and currents 
induced in aircraft systems by lightning are pro- 
portional to the rates of change of the lightning 
electromagnetic fields. Information concerning the 
probabilities of peak rates (changes in current, 
electromagnetic flux density, and other quantities) 


exceeding various limits is necessary to provide the 
basis for determining the susceptibility to lightning 
of electronic systems aboard advanced aircraft and to 
guide the design of protection measures for these sys- 
tems. Therefore, early attention was directed toward 
estimating percentiles of the underlying distributions 
of peak rates based on acquired F-106 time domain 
waveforms (ref. 30). Since then new recording instru- 
ments, called peak detectors, have been added to the 
experimental setup to increase the range of sensitiv- 
ity and to provide new information concerning worst- 
case strikes. Measurements from the peak detector 
give information in the form fa, W{) for i = 1, 2, ..., n, 
where n is the number of flights and each is the 
largest of the peak rates 2 , . that occurs 

as the result of Z^ — z strikes during the ith flight. 

If peak rates are assumed to be statistically inde- 
pendent and to have a common cumulative distribu- 
tion function (e.g., if(x)), then 

P (W t < x\zi) = [H{x)]*t (i = 1,2, ..., n) (9) 

Although the individual peak rates are not observ- 
able, equation (9) directly links ( 2 , w) with the 
underlying distribution. Upon observing a single pair 
( 2 , w) it is known that 2 observations from H(x) fall 
within the range (0, w) and that one or more of these 
is exactly equal to w. To restate this in a slightly 
different way, the data , ( 22 , W 2 ) , (z n , w n ) 

provide the same information as is given by a series 
of largest order statistics in n independent samples 
of varying sizes. 

The work presented herein addresses the problem 
of estimating H(x) for the case in which the form 
of H(x) is assumed unknown. Standard statistical 
techniques, such as relative frequency estimates, are 
not applicable because the maxima are not identi- 
cally distributed. The classic approach (ref. 2) to the 
problem assumes that all realizations of Z \ are quite 
large and concentrates attention on certain limiting 
forms of the extreme value distributions as a basis for 
inference. The more recent theory of concomitants 
of order statistics (ref. 3) considers estimation when 
the members within independent pairs are linked by 
a linear regression model. To our knowledge, there 
are no published accounts which allow a general form 
of the underlying distribution as well as small num- 
bers of observations from the underlying distribution 
which contribute to the maxima. 

In a subsequent section we present maximum 
likelihood estimates of the probabilities that peak 
rates occur within certain intervals. These esti- 
mates are derived from a likelihood function which 
is based on equation (9) and the data (z tl w z ) for 
i — 1, 2, ..., n. We refer to such estimates as 
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nonparametric maximum likelihood estimates be- 
cause they can be calculated without assuming a 
particular form of H(x). 

In addition to providing a nonparametric approach 
to the analysis of data from the peak detector, the es- 
timates might be used to judge the adequacy of vari- 
ous parametric models. They also permit calculation 
of estimates of the probability an aircraft encounters 
strikes of different levels of severity when exposed to 
the lightning environment for various time intervals. 

A Grouped Data Likelihood Function 

The estimation method to be considered assumes 
{Zi,Wj) for i — 1 , 2 ,..., ft are statistically indepen- 
dent pairs and their link with H(x) is as described 
in equation (9). Because P (W^ — 0| Z{ = 0) = 1, the 
joint density of {Zi,Wi) , (Z 2 , W 2 ) , [Z n , W n ) is 

Y[zi{H{wi)] z i h{w { ) n Pi(zi) (10) 
S 1=1 

where S = {i : > 1} and Pi{z) for i = 1,2, ...,n are 

probability mass functions of Z 2 , Z n . Since 
our main interest is H(x) and the second factor 
YlPi{ z i) does not depend on H(x ), we can limit our 

attention to the first factor FI Zi[H{wi)] z i h(wi), 
which is the conditional density of W\, W2 , ..., W n 
given z 1 ,Z2,...,z n . 

The main reason for conditioning is to perform 
an analysis without introducing unnecessary assump- 
tions concerning the form of the densities Pi(z) for 
i — 1,2, If a maximum likelihood estimator of 
H(x) were calculated from the full density in equa- 
tion (10), it would not differ from the maximum 
likelihood estimator if calculated only from the first 
factor. Thus the actual benefit of conditioning is 
that the sampling distribution of the maximum likeli- 
hood estimator depends on the observed z%, Z2, z n 
rather than on the unknown nuisance parameters of 
the densities P{{z) for i = 1,2, ...,n. 

To estimate H(x), consider grouping the maxima 
in a series of intervals having endpoints 0 — ao < 
a\ < ... < ak < dfr+i = 00, all intervals being of 
equal length except for (a*., 00). Let Oj = H(dj) — 
H(aj_i) be the probability of peak rates occurring 
within (a,j_i,aj) for j = 1,2 Observations that 
exceed the full-scale limit of the peak detector are 
easily accounted for by choosing a ^ equal to this 
upper limit. The probability of peak rates exceeding 
a k is 1 ~ 1 — (#1 + #2 + + Ok)' 

Let v 2 = (uii,t; 2 ti---,V(j b+1 ) i ) be a vector of 
indicator variables which define the interval in which 


W{ happens to fall, so that = 1 if a J _ 1 < < 

dj and Vji = 0 otherwise for j = 1,2, ..., k + 1 and i = 
1,2 , n. The joint density of V given 2 is 

/i(v \z,e) = (e z 1 ) vu [(0 1 +8 2 y-e z 1 ] v *' ... 

[(©1 + &2 + ••• + s k) Z ~ (®1 + <*2 + ••• + 8 k - l ) Z ] k ' 

[1 - (0i + 0 2 + ... + 0 k ) z } V(k+1)i (11) 

This density assigns mass to the k + 1 dimensional 
unit vectors (1,0,.. .,0), (0,1,0,...,0), ..., (0,0,. ..,0,1). 

The conditional density of Vi, V 2 , V n , given 
zi,z 2 ,...,z n , is 

<l>=f[fi(.Vi\zi,0 ) (12) 

2= 1 

The grouped data maximum likelihood estimator of 
0 = [0 \ , 62 , Ok) is the statistic 0 — (# 1 , ^ 2 ? • ••> Ok) 
which maximizes equation (12). Calculation of 0 is 
discussed in the next section. 

Maximum Likelihood Estimation 

Consider reparameterization of equation (12) in 
terms of tp = (^ 1 , ^2> • ^ k ) defined by 

6>i +0 2 + ... + Oj = exp U = 1,2 

\ ( 13 ) 

Substituting in equation (12) gives 

n k j 

1 " <t> = 12 ( _2 ») 12 ^3 12 V Qi 

2=1 j=i < 7=1 

n k 

%+l)M [l - ex P (~ z i^j)] (I 4 ) 

1=1 j— l 

The maximum likelihood estimator of ^ is the so- 
lution to dhup/dipj = 0 for j = 1, 2, ..., k. The 
estimating equations are 

n j n 

12 Zt 22 v 9‘ = ^2 z i v {j+l)i [exp( 2 iV>j) - l]~ 

i = 1 <7=1 i = 1 

0=1,2,...,*) (15) 

The main advantage of the reparameterization 
defined in equation (13) is that the ipj can be cal- 
culated by solution of equations in equation (15) 
individually. The solution t/> together with the trans- 
formation in equation (13) gives the maximum like- 
lihood estimator of 0 = (#i, 62, The solutions 
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to equation (15) give estimates of the conditional 
probabilities 


exp(— ipi) = H(ai)/H(a l+l ) ( l = 1,2 


which can be combined to give the product estimator 


H ( aj ) — exp 


k 

-£ 

1=3 


M (j= l,2,...,fc) (16) 


This estimate exists and is unique, providing each of 
the grouping intervals (a^_i, a^) for i = 1, 2 , k + 1 
contains at least one w This is because the expres- 
sions in equation (15) are then positive and finite 
for > 0 and are decreasing functions of xjjj for 
j = 1, 2, k. Under this same condition, 


— <9 2 ln <f>/d\p r dips \ 


= 0 (r # s) 
n 

= 12 z j v (r + l)i expfoVv) 

1=1 

X [exp (Zi^r) - l ]” 2 (r = s = 1 , 2 , 


k) 

(17) 


defines a A; by k diagonal matrix which is positive 
definite; thus, <t> is concave in ^ and the solution to 
equation (15) maximizes (j). 

Asymptotic Properties 

For the independent and nonidentically distrib- 
uted case, reference 31 gives general conditions under 
which a maximum likelihood estimator is consistent 
and asymptotically normal. For the estimators of 
the previous section, the conditions I(i), I(ii), II(ii), 
Il(iii), and III of reference 31 are easy to check under 
the following simpler assumptions: 

A. There exists an upper limit, for example, c, so 

that 1 < < c for i = 1, 2, ... with probability 

equal to 1. 

B. Each 6j for j = 1,2, ...,fc + 1 lies interior 
to (0,1). 

C. JrsW = lim n -1 X -Eld^nfi/d^rd^s) 

i—1 

exists as n — > oo and J = ( J r s{^P)) is positive 
definite. The use of E( ) refers to an expected 
value obtained with reference to the condi- 
tional distribution in equation (12). 

Condition A probably holds because the number 
of strikes that can be recorded over any bounded 
time interval is physically limited. Condition B 
does not differ from what is required for asymptotic 
normality in the standard situation of grouping the 
observations of a random sample. 

Although condition C can hold if Z 2 , Z n 
are not identically distributed, this is the only case 


we attempt to discuss. The expressions on the right- 
hand side of condition C define a diagonal matrix, 
D = ( D rs {'ll))) where 


0 (r^s) 


DrsW { = n 1 £ z i ex P ~ z i S ^ l 

i— 1 y l=r + 1 j 

X [exp [z^r) ~~ 1] _1 (r = s) 


(18) 


The diagonal elements can be written 

n 

n~ 1 J2g r (Z i ) 
i— 1 


where 

g r (z) = z 2 ex p l-z ^ ipi j [exp(z^ r ) - l]” 1 ^ > 1) 

V 1 J 

(19) 

n 

From the law of large numbers, n -1 9r{^i) for 

i=l 

r = 1,2, ...,fc each converge (with probability equal 
to 1) to their expectations which exist (by condi- 
tion A) and are positive (by condition B). This proves 
that condition C holds if Z\, Z 2 , Z n are identically 
distributed. 

Theorem II(iv) of reference 31 states that 
n 1 / 2 (^ — ip 0 ) has a limiting (n — >• 00 ) normal dis- 
tribution with mean vector 0 and covariance ma- 
trix D -1 where D — (D rs (^ 0 )) is evaluated at the 
true parameter vector A standard application 
of the delta method theorem (ref. 32) gives that 
n 1 / 2 (0— 0 O ) also has an approximate normal distribu- 
tion with covariance matrix A -1 , where A= BDB 7 ". 
The elements of B are B rs (0 ) = d'lps/dOr for r and 
s = 1,2, ..., A;, as determined from equation (13). 

Analysis of F-106 Peak Detector Data 

Tables 3 and 4 summarize the largest peak read- 
ings of electromagnetic flux density D-dot and cur- 
rent /-dot as obtained during the 1983 to 1985 flights 
of the F-106 aircraft. On 6 of the first 20 flights, 
maxima of peak D-dot measurements exceeded the 
full-scale limit (48.7 A/m 2 ) of the peak detector. 
This limit was increased to 97.1 A/m 2 and was then 
exceeded only three times in the remaining flights. 
Cases which exceed the full-scale limits are noted 
in table 3. The number of strikes per flight varies 
from 1 to 72, and the typical number is less than 10. 
The largest peak I - dot rate observed in 60 flights 
is 380 kA/fis. 

One way to analyze these data is to consider 
estimates of P(W < x\ z) for z = 1,2, ... obtained by 
grouping the maxima generated by common numbers 
of strikes. This permits weaker assumptions than 
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equation (3), but it greatly expands the dimension of 
the parameter space in that one must then estimate a 
different distribution for each row of tables 3 and 4. 
Further, there are too few observations over much 
of the range of z to give meaningful estimates of 
P(W < x\z) in this way. 

Following the method outlined in previous sec- 
tions, maximum likelihood estimates of 6j = H{aj ) — 
H(cij- 1 ) for j = 1,2 , ...,fc were obtained and are 
given in tables 5 and 6. Standard deviations of the 
estimators are also given in these tables. 

The plots of these estimates in figure 22 are not 
histograms in the usual sense because it is not pos- 
sible to determine the actual percentages of peak 
readings which have occurred in various ranges. 
These plots suggest highly skewed distributions with 
slowly tapering tail areas. 

Because of this feature, moderate increases in 
threshold specifications over much of the upper 
tail range ways by which electromag- 
off in terms of decreasing the probabilities of ex- 
ceeding such specifications. The asymptotic re- 
sults of the previous section give 3.7 ± 1.0 as 
approximate 95-percent confidence limits for the per- 
centage of peak D - dot rates which exceed 48 A/m 2 . 
Similarly, approximate 95-percent confidence limits 
for the percentage of peak /-dot rates which ex- 
ceed 250 kA/fis are 0.7 ±0.2. Actual confidence lev- 
els might differ from 95 percent, depending upon the 
rate of convergence (n — > oo) of the estimators to the 
limiting normal distribution. 

Impact on Airworthiness Criteria 

Although the F-106 research program is of much 
scientific interest, the motivation and objective of 
the research relate to identification and definition of 
lightning hazards to aircraft. These hazards are gen- 
erally thought of as being of two types. The first, 
direct effects, refers to physical damage such as pit- 
ting, burning, magnetic deformation, and destruction 
of dielectric materials such as radomes. The second, 
indirect effects, refers to the hazards caused by elec- 
tromagnetic coupling such as transient voltages and 
currents induced on antennas and cables connecting 
various electronic boxes. Transients can cause dam- 
age or upset of flight- and mission-critical electronic 
systems. Because the F-106 data are primarily use- 
ful for indirect effects, the scope of the remainder of 
this discussion is limited to this topic. In order to 
understand the impact of the lightning research pro- 
gram on airworthiness criteria, it is first necessary to 
discuss how energy couples into electronic systems. 


Internal Coupling of Lightning-Induced 
Electromagnetic Fields 

Electromagnetic coupling is best discussed in two 
parts: external and internal coupling. External cou- 
pling refers to the generation of surface currents 
and charges (or tangential magnetic fields and nor- 
mal electric fields) by a lightning event. Internal 
coupling refers to the generation of electromagnetic 
fields, currents, and voltages on elements of the air- 
craft interior. 

Conceptually, the two parts can be solved in- 
dependently, with the external coupling solution 
defining the sources which then drive the internal 
response. This gives a good approximation if the in- 
ternal coupling does not significantly “feed back” and 
thus alter the external coupling results. This ex- 
ception occurs in instances involving large apertures 
such as an open weapons bay or wheel wells, large 
windows, and nonconducting panels. 

There are three basic 

netic energy penetrates an aircraft interior: aper- 
tures, exposed conductors, and diffusion. Apertures 
refer to openings in the aircraft skin and include such 
items as the cockpit, wheel wells, engine exhaust, 
seams on doors and panels, and areas covered by di- 
electrics. Exposed conductors include antennas of 
various kinds, airspeed probes, and some electrical 
or control cables. Exposed items may receive excita- 
tion either by lightning directly attaching to them 
or from induced fields when lightning is attached 
elsewhere on the aircraft. Diffusion is usually im- 
portant only for carbon fiber composite (CFC) ma- 
terials. Diffusion describes the process by which a 
surface current J s (w) induces an internal tangential 
electric field E s (u) by means of the surface transfer 
impedance Zt{<jj) according to 

E s (uJ ) = Zt(w)J 8 {w) (20) 

where 


sinh k(oj)d 

and uj is the angular frequency, d is the material 
thickness, and and k(u) are respectively the 

intrinsic impedance and the propagation constant of 
the surface material (e.g., CFC). At “late” times (or, 
alternatively, low frequencies), equation (20) reduces 
to the resistive voltage drop along the surface. 

It is useful to summarize the parameters which are 
important to lightning coupling through these pene- 
trations. In order to do so, a simple model of a cable 
is employed and is shown in figure 23. The cable is 
approximated by a transmission line of characteristic 
impedance Zq and having termination loads on each 
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end of impedances Z\ and Z 2 as shown. The electro- 
magnetic coupling sources are approximated in this 
illustration as point sources which include a series 
voltage source Vs and a shunt current source I s as 
shown. In the figure, the transmission line is indi- 
cated schematically by its equivalent circuit model. 

The figure includes a table which lists the pene- 
trations of interest and the lightning coupling vari- 
able upon which the source V s or I s linearly depends. 
For example, the source V s for an open hole aperture 
is a constant time B-dot. Contained in parentheses 
is the source variable which relates to the lightning 
environment (/-dot in this example). 

It should also be pointed out that the table does 
not, by itself, give the entire picture with respect 
to what is important for coupling. For example, 
if Z\ and Zj are small or zero and the cable is 
short with respect to wavelengths in the lightning 
pulse (which is often the case), then for aperture 
coupling dominated by the B field the current I\ is 
proportional to B (or /, because /1 for this situation 
is simply proportional to the time integral of V$). 
Therefore, even though V s is proportional to /-dot, 
I\ is proportional to /. On the other hand, if Z\ and 
Z 2 are large, then V\ relates to /-dot and not to /. 

In addition, in some cases where the transient 
responses of the cable are related to /, then the 
total waveform of / is important. This also requires 
quantification of the energy delivered to a load, which 
places a requirement on the action integral (ref. 33) 

f°° 1 2 dt 
Jo 

This would be particularly important for damage 
of electronics as well as upset. Therefore, even for 
aperture coupling, /-dot, /,D- dot, and the action 
integral are important. For diffusion, only I and the 
action integral are important. 

Certain other knowledge of lightning environ- 
ments is required for system considerations. For ex- 
ample, upset of a digital system is determined by the 
amplitude of a pulse, the pulse width, and the repeti- 
tion rate of pulses. Therefore, specifications on these 
parameters are also of interest. 

An example of a comparison of a calculation 
and a measurement on an internal wire based on 
the model of figure 23 is shown in figure 24. The 
wire is inside the F-106 fuselage and is 9 m long, 
shorted to the structure on one end and terminated in 
a 50-fi resistance on the other end. The wire passes 
within 0.61 m of some closed wheel well doors whose 
seams were modeled in the analysis as inductive 
seams because of their observable large gaps. The 
model for calculating the wire response is identical to 
that given in figure 23. The gap transfer inductance 


used in the model for the wheel well door gaps was 
chosen to be 6.45 nH-m. The lightning nose current 
I is given in figure 25. By inspection, it is clear 
that the large oscillations dominating the current 
response are related to the derivative of /. This 
is reasonable, since the wire inductance is 1.5 /iH 
and the L/R time constant with the 50-fi load is 
only 30 ns, which is much shorter than the period 
of oscillation (ss250 ns) of the cable current. In this 
frequency regime, the current response is limited by 
the 50-fi load impedance, which is a large impedance 
compared with the cable inductance. Therefore, 
the current (or voltage) would be proportional to /, 
which it clearly is. 

Note that this is an example of an internal wire 
response caused by a lightning event of low peak cur- 
rent and low /-dot (^20 kA/^s peak). Even so, it 
produces an internal current pulse having frequencies 
related to aircraft and cable resonances. Aircraft ex- 
ternal transient responses due to lightning have spec- 
tra with resonance peaks nominally in the range of 1 
to 50 MHz, which will excite internal responses close 
to clock frequencies of aircraft computers (ref. 8). 

Impact on Testing 

The main issues of interest in testing include /- 
dot, D- dot, the temporal durations of /- dot and 
D- dot, and the late time currents, which are im- 
portant for diffusion through CFC structures. By 
“late” time current we mean times exceeding a diffu- 
sion time constant. For times shorter than this time 
constant, a composite behaves similarly to a metal. 
For later times, current redistribution and diffusion 
through the composite structure is observed. Ap- 
proaches which have been used to increase both /- 
dot and D-dot capabilities of generators for aircraft 
testing include the use of discrete peaking capaci- 
tors (refs. 34 to 36), distributed peaking capacitors 
(ref. 37), high-voltage oil-filled low-inductance Marx 
generators (refs. 38 and 39), and the so-called shock 
excitation techniques (ref. 40). In order to increase 
late time currents, a crowbar switch is used to pro- 
vide for an exponential decay of a Marx generator 
as opposed to the normal damped sinusoid (ref. 41). 
It is clear that there is no consensus yet on the best 
way to accomplish these types of tests. 

The maximum F-106 peak recordings of current 
and current rate of change are shown in table 7 
along with a summary of the recent past criteria 
and new criteria established in the joint Atmospheric 
Electricity Hazard Protection Program for Aircraft 
(AEHP) (ref. 42). (The AEHP is a multiagency 
program managed by Wright Aeronautical Laborato- 
ries to develop and demonstrate lightning protection 
technology for digital systems onboard composite 
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aircraft.) Note that the absolute largest F-106 cur- 
rent rise rate measurement of 380 kA / fxs is consider- 
ably greater than the old or new criteria for current 
rate of change. The approximate 95-percent confi- 
dence limits for the percentage of peak /- dot rates 
which exceed 250 kA///s are 0,7 ± 0.2, as shown in 
figure 22. Considering that criteria are usually es- 
tablished at about the 99th quantile, 200 kA//is as 
used in the AEHP program is appropriate. 

The impact of the new criteria on testing for 
/-dot, D-dot, and their temporal durations is not 
fully known, but some general observations can be 
made. In order to increase I - dot or D-dot and 
to lengthen their duration, either peaking capaci- 
tors with high-inductance Marx generators or with 
low-inductance oil-filled Marx generators should be 
used. A low- voltage capacitor bank can possibly also 
be used but the transmission line return conductor 
system must be carefully designed and must be very 
close to the aircraft under test (ref. 43). The relative 
expense of providing large values of these parame- 
ters increases greatly as the size of the test object 
increases. For example, testing a vehicle the size of 
a Boeing 747 at /-dot = 200 kA//zs is not within the 
state of the art, but this type of test can be done on 
a small fighter-size aircraft. 

Concluding Remarks 

This research effort has resulted in the first sta- 
tistical quantification of the electromagnetic threat 
to aircraft based on in situ measurements. Previous 
estimates of the in-flight lightning hazard to aircraft 
were inferred from ground-based measurements. The 
electromagnetic measurements made on the F-106 
aircraft during these strikes have established a statis- 
tical basis for determination of quantiles and “worst- 
case” amplitudes of electromagnetic parameters of 
rate of change of current and the rate of change of 
electric flux density. The 99.3 percentile of the peak 
rate of change of current on the F-106 aircraft struck 
by lightning is about two and a half times that of 
previously accepted airworthiness criteria. The find- 
ings are at present being included in new criteria con- 
cerning protection of aircraft electrical and electronic 
systems against lightning. Since there are at present 
no criteria on the rate of change of electric flux den- 
sity, the new data can be used as the basis for new 
criteria on the electric characteristics of lightning- 
aircraft electrodynamics. In addition to there be- 
ing no criteria on the rate of change of electric flux 
density, there are also no criteria on the temporal 
durations of this rate of change or rate of change of 


electric current exceeding a prescribed value. Results 
on pulse characteristics presented herein can provide 
the basis for this development. The newly proposed 
lightning criteria and standards are the first which re- 
flect actual aircraft responses to lightning measured 
at flight altitudes. 

The focus of the discussion on phenomenology 
and modeling of lightning strikes centers on triggered 
lightning, since most of the strikes to the F-106 are 
triggered events. In every case of 49 ultrahigh fre- 
quency radar observations of the F-106 being struck 
by lightning, the data indicate that the channel be- 
gins at the aircraft and travels away from it. This 
evidence prompted the development and application 
of nonlinear air breakdown modeling of the aircraft 
triggering which gives results in good agreement with 
measured data. (Subsequently it was found that 
most, but not all, of the F-106 strikes were triggered.) 
Given information on lightning attachment location, 
the injected current waveform, etc., the linear mod- 
eling technique developed provides results for sensor 
responses over the surface of the aircraft which are 
in excellent agreement with the measured data. This 
technique is appropriate as a first step in the com- 
plex problem of prediction of internal induced effects 
in arbitrary aircraft types. 

The discussion on maximum likelihood estimation 
of the distribution of peak rates of current and elec- 
tromagnetic flux density on the basis of the F-106 
peak detector data is the first known attempt to es- 
timate a general form of underlying distribution on 
the basis of observing only the maxima of random 
variables. The main thing that this analysis provides 
is estimates of the probabilities that peak rates oc- 
cur in various ranges. These estimates suggest highly 
skewed distributions with slowly tapering tail areas, 
which are of much interest in evaluating the lightning 
threat to aircraft. 

Finally, the research reported herein has con- 
tributed to basic scientific knowledge by establishing 
new methods and results quantifying the electromag- 
netic interaction of lightning with aircraft. It must 
be recognized, however, that these results were ob- 
tained at altitudes ranging from 15000 ft (4570 m) to 
about 40000 ft (12 190 m), where strikes are quite ob- 
tainable. Since lower altitude strikes are much more 
difficult to obtain, much work remains to quantify 
low-altitude lightning-aircraft electrodynamics. 

NASA Langley Research Center 
Hampton, Virginia 23665-5225 
April 27, 1987 
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Appendix A 

Air Chemistry Coefficient Formulas 
(Ref. 44) 

Calculation of E re \: 


^rel — 


•"mag 


Pr (1 + 2.457P 0 - 834 ) 


E, 


< 0.07853 (l + 2.457P 0834 ) j (Al) 


E 

E t 


E Te] = EEM. (-1.195P 0 834 ) 


^ > 3.015 + 1.195P a834 ) (A2) 

Pr J 


^rel — 



E m ag - / 0.6884P 0 ' 834 ' 

)1 

1/2 ^ 
0.688P 0 - 834 ( 


. * + V 2 

/ 

2 j 

(allother ^ raag ) 


0 


where P is the percent water vapor, p r is the relative 
air density, and P m ag is in electrostatic units (ESU). 

Calculation of electron attachment rate a e : 

100- P 


<*e 


100 


[0:3(1 + 0.344P) + <* 2 ] 


where 


Q 2 = 1.22 X 10 8 p r exp ( — 21.15/Prel) 


a 3 = p 2 (6.2 x 10 7 + 8.0 x 10 10 P 2 el ) / 
1 + lO 3 ^] [pfel ( 1 + 0.03E 2 e| j 


Calculation of avalanche rate G: 

G = 5.7 x 10 s p r Y b /(l + 0.3 Y 2 - 5 ) 
where Y = E re \/ 100. 

Calculation of electron-ion recombination coefficient 
/? and ion-ion neutralization coefficient 6: 

(3 = 2 x 10 -13 + 2.8 x 10 _12 (P) 1/3 m 3 /s 
6 = 2 x 10“ 13 + pr 2.1 x 10 _12 m 3 /s 

Calculation of electron mobility p € : 

lOOpta 


Me = 


100 - P + PR 


where 


R = 1.55 + 210/(1 + 1 1.8£ re i + 7.2£ 2 el )(m/s)/(V/m) 
and 

p a = {[(16.8 + P rel )/(0.63 + 26.7P rel )] 0 ' 6 } /3p r 


Calculation of ion mobility pf 
2.5 X 10“ 4 


Pi = 


Pr 


(m/s)/(V/m) 
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Appendix B 

Aircraft Enhancement Factors 

One of the uses of modeling is to determine the 
enhancement factors of the metallic F-106 aircraft 
if it is placed in a static electric field or has a net 
charge on its surface. These enhancement factors 
can be defined as the ratio of the local electric field 
at any location on the aircraft to the ambient field 
magnitude in the absence of the aircraft. Of course, 
for the case of a net charge on the aircraft, the 
enhancement factor must be defined as the ratio of 
the local electric field to the net charge. 

The importance of enhancement factors is that 
they determine the conditions under which triggered 
lightning can occur. The local electric field at any 
point on or near the aircraft must exceed the value 
necessary for air breakdown in order for this to 
happen. Therefore, the enhancement factors can be 


used to find the minimum ambient field for which 
triggering will occur at a given altitude. 

The calculated enhancement factors at selected 
locations on the block model of the F-106 are shown 
in figure 26 for three directions of ambient electric 
field and a net positive charge on the plane. The 
fields in figure 26 are in volts per meter for an 
ambient field of 1 V/m or in volts per meter for a 
1-//C net charge. 

The calculated fields are not the largest fields that 
would be seen on the actual F-106. Because of the 
finite spatial resolution of the model, details such as 
the nose boom are impossible to resolve. Hence, the 
large fields around the tip of structures such as this 
are not seen in the model. However, these large fields 
exist only in a very small volume, and whether or not 
they are the controlling factor in triggering lightning 
is unknown. A large volume of smaller fields with a 
higher total energy may be necessary for triggering 
to occur. 
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Table 1. Comparison of High- and Low-Altitude F-106 Triggering Fields as a Function of Charge and Orientation 


Aircraft 

charge 

Magnitude* 1 , kV/m, for electric field orientation of — 

Back to 
front 

Front to 
back 

Right to 
left 

Left to 
right 

Bottom to 
top 

Top to 
bottom 

0 

190 

190 

280 

280 

470 

470 


(250) 

(250) 

(380) 

(380) 

(590) 

(590) 


107 

200 

200 

200 

300 

600 


(140) 

(230) 

(210) 

(210) 

(350) 

(660) 

—0.5 Qm 

273 

150 

200 

160 

640 

300 


(230) 

(160) 

(210) 

(210) 

(660) 

(350) 

Qrn 

24 

50 

20 

20 

80 

80 


(34) 

(60) 

(40) 

(40) 

(120) 

(160) 

~~Qrn 

24 

15 

20 

20 

80 

80 


(60) 

(34) 

(40) 

(40) 

(160) 

(120) 


a Low-altitude field magnitude is in parentheses. 




Table 2. Trigger Conditions for Aircraft of Various Sizes and Shapes 


Aircraft 

Field 

orientation 

Qm, 

mC 

■E'mag) 

kV/m 

Half-size F-106 

Nose to tail 

0.45 

100 


Top to bottom 

.45 

250 

Double-size F-106 

Nose to tail 

7.16 

100 


Top to bottom 

7.16 

250 

Normal-size F-106 with 

Nose to tail 

2.00 

120 

straight wing 

Top to bottom 

2.00 

240 

Normal-size C-130 

Nose to tail 

6.90 

190 


Top to bottom 

6.90 

210 

Normal-size F-106 

Nose to tail 

1.79 

100 


Top to bottom 

1.79 

250 


2 



Table 3. Largest Peak Rates of Change of Electric Flux Density for 62 Flights of the F-106 Aircraft 


Number of 
strikes 


Largest peak rates, A/m 2 


4.5, 9.9, 15, 15, 17, 17, 17, 18, 20, 27, 43, 47, a 48 
7, 9.7, 15, 23, 24, 26, 29, 33 , 39 

12, 24, 26, 28, 28, 31, 40, a 48 

13, 39 

17, 26, 28, 28, 28, 40, 75 
16, 17, 37, 40, 78, a 97 
30, 37, 42, 42, a 97 
27, 30, “48 


“Certain peak rates exceeded the full-scale limit of the peak rate detector. Initially 
this limit was 48.7 A/m 2 and was later increased to 97.1 A/m 2 . 









Table 4. Largest Peak Rates of Change of Current for 60 Flights of the F-106 Aircraft 


Number of 
strikes 

Largest peak rates, kA///s 

1 

1.90, 7.6, 9.6, 10, 20, 20, 42, 54, 78, 82, 96, 134, 136 

2 

7.6, 11.4, 24, 28, 34, 40, 84, 88 

3 

1.91, 1.91, 20, 78, 88, 300 

4 

7.6, 42 

5 

5.8, 11.4, 44, 48, 60, 68, 198 

6 

3.8, 15.2, 22, 50, 110, 162 

7 

5.8, 22, 30, 76, 98 

8 

86, 134, 240 

9 

194. 380 




Table 5. Maximum Likelihood Estimates of the Probabilities That 
Peak Rates a of Electric Flux Density Occur Within Selected Ranges 


Lower grouping 
limit, A/m 2 

Number of 
flights 

Number of 
strikes 

Estimate 

Standard 

deviation 

0 

i 

i 

0.098 

0.057 

6 

3 

5 

.197 

.068 

12 

10 

27 

.314 

.063 

18 

4 

8 

.070 

.016 

24 

13 

47 

.162 

.027 

30 

6 

43 

.048 

.008 

36 

9 

66 

.051 

.007 

42 

5 

25 

.023 

.003 

48 

11 

177 

.037 

.005 


a For 
and are 


example, maxima of peak rates are counted in the range of 6 to 12 if they are less than 12 A/m 2 
greater than or equal to 6 A/m 2 . 
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Table 6. Maximum Likelihood Estimates of the Probabilities 
That Peak Rates“ of Current Occur Within Selected Ranges 


Lower limit, 
kA/fis 

Number of 
flights 

Number of 
strikes 

Estimate 

Standard 

deviation 

0 

30 

135 

0.876 

0.012 

50 

15 

133 

.079 

.011 

100 

5 

21 

.017 

.002 

150 

4 

47 

.013 

.002 

200 

3 

50 

.008 

.001 

250 

3 

32 

.007 

.001 


“For example, maxima of peak rates are counted in the range 50 to 100 if they are less than 100 kA///s 
and are greater than or equal to 50 kA///s. 
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Table 7. Comparison of F-106 Data With Existing Criteria 
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Figure 1. Location of electromagnetic sensors. 












;ure 5. Rate of change of electric flux density from sensor located under forward fuselage. (Recorded at 
10-ns sample interval.) Waveform recorded simultaneously with that in figure 6. 
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Figure 7. Model of F-106 used in finite-difference computer code. 






















Figure 11. Comparison of nonlinear model calculated response (bottom to top, jEVnag = 0.64 MV/m, 
Q = Q ml 2) with measured data for B-dot longitudinal sensor. 
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Edge of problem space 



I 

(a) Natural lightning. 



(b) Triggered lightning. 

Figure 12. Model showing locations of current source used in computer code. 
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Figure 21. Dimensional d 
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(c) Normal-size F-106 with conventional straight wings. 
Figure 21. Continued. 


I 
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(d) Normal-size C-130. 
Figure 21. Continued. 
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ELECTRIC FLUX DENSITY , A/m k A / p s 


Cells not illuminated by A cell illuminated by lightning- Cells not illuminated by 

lightning-induced fields induced fields coupling through lightning-induced fields 

a penetration 



Lumped parameter model of a transmission line with characteristic impedance 
Z Q . Each cell exposed to lightning-induced fields has a voltage source V s and a 
current source l s . The sources are dependent upon the lightning environment 
as shown below for various penetrations. 


Source Relationship to Lightning Environment 


Penetration Tvpe 

V s Proportional to 

L Proportional to 

Open hole aperture 

B-dot (l-dot) 

D-dot (1) 

Resistive seam aperture 

B (1) 

negligible 

Inductive seam aperture 

B-dot (l-dot) 

negligible 

Diffusion through CFC 

Z t B (1) 

negligible 

Exposed conductor - 

Direct attachment 

negligible 

1 

Exposed conductor - 
Lightning strikes 

B-dot (l-dot) 

D-dot (1) 


aircraft nearby 


Figure 23. Transmission line model of lightning coupling to cable. 






? or internal wire response of figure 24. 
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